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Abstract 

In this paper a greedy algorithm to detect conflict cliques in interval 
graphs and circular-arc graphs is analysed. In a graph, a stable set re- 
quires that at most one vertex is chosen for each edge. It is equivalent 
to requiring that at most one vertex for each maximal clique is chosen. 
We show that this algorithm finds all maximal cliques for interval graphs, 
i.e. it can compute the convex hull of the stable set polytope. In case 
of circular-arc graphs, the algorithm is not able to detect all maximal 
cliques, yet remaining correct. This problem occurs in the context of rail- 
way scheduling. A train requests the allocation of a railway infrastructure 
resource for a specific time interval. As one is looking for conflict-free train 
schedules, the used resource allocation intervals in a schedule must not 
overlap. The conflict-free choices of used intervals for each resource cor- 
respond to stable sets in the interval graph associated to the allocation 
time intervals. 

1 Introduction 



1.1 Motivation 

One important aspect of railway optimization is scheduling trains in a main 
station area [6, 2]. In such an area the track topology is complex, and many 
different routes to travel between the area entrances (called portals) and plat- 
forms exist. Here, an appropriate assignment of exactly one of these routes to 
each train is crucial. The choice of these routes has to be conflict-free, i.e. each 
infrastructure resource can be occupied by at most one train at the same time. 
Hence, the allocation time intervals at each resource belonging to the chosen 
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train routes can not overlap ([9]). These restrictions over the concurrent al- 
location time intervals at each resource transforms into finding stable sets in 
corresponding interval graphs. The stable set problem itself can be modeled by 
formulating an integer linear program (ILP) and solving it using a commercial 
solver, once an objective has been formulated. Since a good ILP formulation 
has a strong relaxation which speeds up the solution process, we look for an 
efficient ILP description of the stable set problem in an interval graph. 

1.2 The basic idea 

Let consider the track topology, consisting in a set of resources, as described 
in [7]. Each route/time assignment for a train allocates the used resources 
over certain time intervals, depending on signal positions, train dynamics, track 
topology, and additional safety regulations. A feasible train schedule assigns a 
route and times to each train such that no resource is allocated at the same 
time, thereby guaranteeing a conflict-free schedule. 

Several approaches exist in the literature for modeling conflicting allocations 
over resources. The conflict modeling in the conflict graph approach of [10] is 
simple: each available route/time assignment of a train corresponds to a node 
in the conflict graph. Each time two nodes of distinct trains would allocate a 
resource at the same time, a conflict edge between these nodes is introduced 
into the graph. Additionally, the fact that a train is only allowed to have one 
route/time assignment is modeled by interconnecting all nodes of the same train 
by edges thereby forming a clique. A solution of the train routing problem in 
this model is a stable set with cardinality equal to the number of trains n = \T\, 
where each chosen vertex assigns a route/time to the corresponding train. As 
it is a stable set of cardinality n, this assignment is guaranteed to be conflict- 
free and each train gets exactly one route/time assigned. We can illustrate this 
conflict modeling concept by drawing all the allocation time intervals of a single 
resource during a period of time (see later Figure 1). This leads to the following 
ILP formulation for all the stable sets (feasible solutions) of the conflict graph G, 
where the variable x^ correspond to the node ij E G, meaning the route/time 
assignment j for train i: 



m{i) 

= 1 for alH = 1, . . . , n (1) 

i=i 

Xik + Xji < 1 for all r ik w r jt (2) 
Xij €{0,1} (3) 

where r ik Tj l means that assignment k of train i is in conflict with assignment 
I of train j. 
Let 

STABl(G) := {x\x >0,x satisfies (1), (2)} 
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be the polytope of the non-negative vectors x satisfying (1) and (2). It is easy 
to see that integral solutions fulfilling (1), (2), and (3) are exactly the incidence 
vectors of stable sets of nodes of G. 
Let denote with 

STAB(G) := conv{x s G {0, 1} V \S C V is stable set} 

the convex hull of the incidence vectors of all stable sets of nodes of G. Clearly, 
STAB(G) C STAB1(G) and the integer points contained in both sets are the 
same. The two polytopes are equal if and only if the graph G is bipartite and 
has no isolated nodes [5] . This is in general not the case when looking at train 
scheduling problems. Hence, we have to improve on the STAB1(G) formulation. 

1.3 Gathering conflicts 

In the following we restrict our view to a single resource of the track topol- 
ogy, consisting in a set of topology elements [7] . Instead of looking at pairs of 
overlapping allocation time intervals like in the conflict graph approach, [4, 1] 
introduces a more sophisticated attempt which gather them into groups of con- 
flicting intervals. All possible assignments using the resource at the same time 
are grouped in a conflict clique, and at most one of these can use the resource 
(see later Figure 2). 

Let the graph G = (V, E) denote the induced interval graph from the alloca- 
tion time intervals for this resource, which is built as follows: for each allocation 
time interval a vertex v G V is created and two vertices are connected with an 
edge e G E if the two intervals intersect. On the other hand, in periodic schedul- 
ing we deal with circular time axes and the induced intersection graph G is called 
circular-arc graph. The set of allocation time intervals is called arc model. 

A feasible solution of the train scheduling problem should be conflict-free in 
each resource. It means that the set of chosen assignments represent a stable set 
in the interval graph of each resource. The conflict graph formulation describes 
the set of all stable sets in the graph G using the incidence vectors and (integer) 
linear constraints avoiding the choice of two adjacent nodes. Ideally, one would 
be able to describe efficiently the convex hull of the incidence vectors of all stable 
sets in G as a polytope. Since the polytope is naturally integer it induces an 
efficiently solvable ILP formulation. Chvatal [3] proved that one can describe 
efficiently this convex hull. This description imposes that the amount of stable 
set vertices over each maximal clique of the interval graph should be at most 
one. 

Section 2 explains the greedy algorithm to detect the conflict cliques and 
in Section 3 we prove that this algorithm finds all maximal cliques for interval 
graphs (i.e. describes the convex hull). 
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2 A greedy algorithm for grouping conflicts into 
cliques 



In this section we explain the algorithm introduced by [4, 1] for detecting con- 
flicts and grouping them into cliques. 

2.1 Allocation Schema 

For an easier understanding of the algorithm we first introduce a graphical 
representation of resource allocations called Allocation Schema (AS). We create 
for each resource an AS, which is a representation of the resource allocation 
intervals by the potential assignments over time. We chart the allocations of 
the different assignments by a horizontal line with the start of the allocation as 
the offset and the allocation duration as the length of the line. An example of 
AS is given in Figure 1. 
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Figure 1: In the Allocation Schema (AS) the allocations of a resource by train 
assignments are represented by horizontal lines. 



2.2 Algorithm for finding conflict cliques 

To find all conflict cliques we suggest an algorithm which operates on the AS 
and not on the intersection graph directly. Algorithm I basically consists of two 
major steps: 

1. Sort the start and end times of the allocation time intervals according to 
time. 

2. Walk through the times, keep track of the currently open intervals in a list 
and construct the set of conflict cliques among the currently open intervals 
at the end times. 

Note that we do not have to form a clique at the end time of each interval, but 
only at end times where a new interval has been opened since the last iteration. 
Additionally, we avoid adding a clique if it just consists of one interval since it 
would not be a necessary constraint. A result of Algorithm 1 is illustrated in 
Figure 2. 



In 
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Algorithm 1 Conflict Cliques 



of a resource 



i}, boolean 



Input: Set of all allocation time intervals {[s(«), /(*)]}. —1 
Output: Set C r of conflict cliques 
1: Create a list L of 2n tupels (time 6 R + , identifier € {1,. 
is-endtime) 

Sort L according to the first key time, for same time starttimes before end- 
times 

Initialize list of currently open intervals := 
Initialize new.starttime := FALSE 
for i = 1 to In do 

if is_endtime = FALSE then 
:= O U {identifier 1 } 
newstarttime := TRUE 



else 



if [newstarttime 
C r := C r U 
newstarttime := 

= \ {identifier } 
return C r 



TRUE)kk{\0\ > 1) then 
= FALSE 



I " 
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Figure 2: Several train itineraries allocate a resource during different time in- 
tervals (horizontal bars). The dashed vertical bars illustrate conflict sets among 
the allocation intervals formed by Algorithm 1. 



2.3 Complexity and remarks 

Algorithm 1 runs in O(rtlogn) where n is the number of allocation intervals. 
In a first step we have to sort the allocation intervals according to time which 
takes O(nlogn) time [8]. We then loop 2n times where the steps inside the loop 
take O(l) time which results in a running time complexity of O(n) for the loop 
and O(nlogn) for the whole algorithm. 

For non-periodic railway scheduling problems the described algorithm is op- 
timal, i.e. it finds all maximal cliques in the corresponding interval graph and 
thereby also the convex hull of the stable set incidence vectors. The proof fol- 
lows in Section 3. In periodic railway scheduling problems, the allocations of a 
resource reoccur after a fixed period length and hence allocation intervals may 
overlap the period length. To deal with these overlaps of the period length, we 
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can use the modulo function to draw overlapping intervals by two horizontal 
lines. The first line starts at time and ends at the original end time of the 
allocation modulo the period length. The second line starts at the original start 
time and ends at the period length. 

There are two drawbacks in this periodic setting: the first is the unknown 
quality of the polytope described by the maximal clique. For instance, in the 
chordless odd cycles K2k+i for k > 2 the non-integer solution where all vertices 
have assigned the value ^ is an extremal point of the polytope described by 
the maximal cliques, making it in general not the convex hull for the incidence 
vectors of a stable set. It is unclear how one can find efficiently the convex 
hull of the stable set incidence vectors in the periodic setting and how good the 
maximal clique description approximates this convex hull. 

The second problem is, that the current greedy algorithm fails to detect all 
maximal conflict cliques. For instance, in the example of Figure 3 the described 
method fails to detect a conflict clique: Typically, in practice the length of 
a period is one hour and the allocation time intervals are not longer than 3 
minutes. Hence these special cases where a conflict clique is not detected by the 
algorithm most likely will not occur. 
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Figure 3: Special case, where the three clique consisting of the intervals A,B, 
and C is not detected 



3 Optimality proof for the non-periodic case 

In this section we present a proof that the greedy algorithm is optimal in case 
of non-periodic scheduling. 

Definition 3.1. For a time interval i of a resource r, denote its starting and 
finishing time by s(i) and f(i), respectively. An interval i is called active at 
time t iff s(i) < t < f(i). 
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Property 3.2. A subset Jc of all allocation intervals V of a resource r forms 
a clique Cj of the associated interval graph G r , iff there exists a time segment 



s(J c ) ■■= max {s(i)} ,f(J c ) ■= mm {f(i)} 



during which all intervals i £ Jc are active. Furthermore, for any t < s(Jc) or 
t > f(Jc), at least one of the intervals i £ Jc is not active. 

Theorem 3.3. Algorithm 1 computes a minimum edge clique cover for the 
interval graph G r associated with the allocation intervals I r . Each clique of the 
edge clique cover is maximal. 

Proof. Clearly, C r computed by Algorithm 1 is an edge clique cover. Further- 
more, every clique C £ C r is maximal. Suppose some C £ C r was not. Then, 
there would exist an interval a Jc, conflicting with every interval b £ Jc- 
But because the clique C, which was added to C r in line 11 of Algorithm 1, 
is equal to the set Oc of all open (i.e., active) intervals at time t — f(Jc), 
and a f Oc, it must hold that either s(a) > f(Jc) or /(a) < s(Jc), where 
s( Jc) — s (identifier*), the starting time of interval identifier 1 defined in line 7 
of the algorithm. By Property 3.2, both cases contradict the assumption that 
a conflicts with every interval b £ Jc- 1 

Finally, to show that I :— \C r \ is minimum, it suffices to identify a subset S 
of the conflicts that requires at least I cliques to be covered. This subset can 
be constructed inductively as follows. Denote by d, 1 < i < I the cliques of 
C r in the order in which they are found by the algorithm. For every clique Ci, 
denote by Si, fi intervals for which s(Ci) = s{si) and f(Ci) = /(/i), respectively. 
Let a be the conflict between Si and fi, if Si ^ fi, or an arbitrary but fixed 
conflict between Si — fi and another interval a £ Ci. In either case, for any 
i < j, conflicts Ci and Cj cannot be covered by the same clique, since f{fi) < 
s(sj), and hence, intervals fi and Sj are not in conflict. Therefore, the set 
S = {c^ 1 1 < i < £} is a subset of conflicts that are pairwisc not covcrable by 
the same clique. Hence, any edge clique cover must consist of at least \S\ = I 
cliques. □ 

Corollary 3.4. Algorithm 1 finds all maximal cliques in the interval graph G r . 

Proof. Suppose there exists a maximal clique M C r . Clearly, M can neither 
be a superset nor a subset of any C £ C . Hence, M can only contain subsets 
of intervals of two consecutive cliques Ci, Ci + i £ C r , i.e., M — C • U C- +1 , with 
nonempty subsets C[ C Ci and C' i+1 C Cj+i. Now, M (fi d can only hold if 
there is an interval a in C' i+1 that is not in C[. M is only a clique, if a is in 
conflict with every b £ C[. But if this is the case, then M C Ci+i, since all 
intervals that are open at s{a) are, by Algorithm 1, in Cj+i. □ 



1 In the algorithm, a suitable tie break should be chosen, so that for intervals i, j with 
f(i) = s(j), the starting time of j is considered before the finishing time of i. 
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